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ABSTRACT 



An 11-day monitoring campaign in late 2005 reveals clear correlation in po- 
larization between the optical emission and the region of the intensity peak (the 
"pseudocore" ) at the upstream end of the jet in 43 GHz Very Long Baseline 
Array images in the highly variable quasar PKS 0420—014. The electric-vector 
position angle (EVPA) of the pseudocore rotated by about 80° in four VLBA ob- 
servations over a period of nine days, matching the trend of the optical EVPA. In 
addition, the 43 GHz EVPAs agree well with the optical values when we correct 
the former for Faraday rotation. Fluctuations in the polarization at both wave- 
bands are consistent with the variable emission arising from a standing conical 
shock wave that compresses magnetically turbulent plasma in the ambient jet. 
The volume of the variable component is the same at both wavebands, although 
only ~ 20% of the total 43 GHz emission arises from this site. The remainder 
of the 43 GHz flux density must originate in a separate region with very low 
polarization. If 0420—014 is a typical case, the nonthermal optical emission from 
blazars originates primarily in and near the pseudocore rather than closer to the 

^ , central engine where the flow collimates and accelerates. 
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Introduction 



Our knowledge of the jets of active galactic nuclei (AGN) becomes less certain as we 
consider regions closer to the central engine. However, these regions are the sites of especially 
interesting physical phenomena, such as acceleration and collimation of the jet flow, forma- 
tion of the observed "core" of the jet, energization of relativistic electrons, and changes in 
the geometry of the magnetic field. While it is widely accepted that synchrotron emission is 
the main mechanism behind continuum production, establishing the location of this emission 
at different wavebands has been problematic. Imaging with very long baseline interferom- 
etry (VLBI) can resolve different components of the jet at radio wavelengths. In addition, 
variability of the flux density and spectral shape provides information on the small-scale 
behavior at all wavebands, but the inability to produce high quality images at submilliarc- 
second resolution at wavelengths shorter than 3 mm allows multiple interpretations that are 
compatible with the data. 

There is, however, an emerging method for establishing the location of the emission at 
short wavelengths. Similarities of unique polarization signatures of features in both VLBI 
images and opti cal measurement s proy ide a possible conne c tion b etween the two wavebands. 
Recent work by Gabuzda et ali ( 2006 ) and Jorstad et al. (12007 ). that ex tend earlier stud- 
i es by IWardle. Moore, fc Angel ( 1984T). Isitko. Schmidt, k Steinl jl985h . iGabuzda fc Sitko 
( 119941 ). and lLister fc Smith! (120001 ). reveals a correlation between VLBI pseudocore (see be- 
low) and optical polarization angle in a number of blazars. These results suggest that the 
two emission regions are either cospatial or sepa rate b u t with both threaded by magnetic 
fields having similar geometries. [Here we follow I Jonesl (119881 ). in referring to the compact 
feature at the upstream end of the jet in VLBI images as the "pseud ocore" . This may be an 



actual stationary feature, such as a standing shock wave (see, e.g., iDaly fc Marscherlll988 



Cawthornd 120061), or the point in the divergin g jet flow where the optical depth is of order 
unity falandford fc Konr3ll979l : lKonig]||l98lh .] 



The most direct way to establish cospatiality of the emission at different wavebands is to 
observe coordinated variability of polarization that is seen both in the VLBI images and in 
measurements at shorter wavelengths. In this Letter we report observations of polarization 
variations in both the optical emission and the pseudocore as measured o n 43 GHz Very Long 
Base line Array (VLBA) images of the quasar PKS 0420-014 (z = 0.915; Hewitt fc Burbidge 
19931 ) that span nine days, with no significant time lag. These results place the nonthermal 
optical emission within the 7 mm pseudocore of the jet. 
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Observations 



We collected optical data at two sites. At the Steward Observatory 1.55 m Ku iper tele- 
scope, we used the SPOL spectropolarimeter ([Schmidt. Stockman, k. Smith! Il992l ). making 
a total of nine measurements from 2005 October 23 to November 3. A 600 lines/mm grating 
covered an entire diffraction order of ~ 400-800 nm at a resolution of 1.7 nm. The spectropo- 
larimetry used slit widths of either 2" or 3" along with a 10" spe ctral extraction aperture . 
Differential spectrophotometry relative to field stars calibrated by ISmith fc Balonekl (119981 ) 
emp loyed a 5.1" x 12" ( or 14") aperture. Details of the reduction procedure can be found in, 
e.g., ISmith et al.l (120031 ). We also performed polarimetry and photometry at R band on the 
AZT-8 70 cm telescope at the Crimean Astrophysical Observatory. The signal-to-noise ratio 
was low, however, owing to the faintness of 0420—014 (R = 17.8-18.3) during the campaign. 
We therefore do not display the data here, but note that they are consistent with the Steward 
Observatory results. 

We observed 0420—014 at 43 GHz with the VLBA at four epochs during the campaign: 
October 24, 28 and November 1, 2. After correlation at the Array Operations Center of the 
National Radio Astron omy Observatory (N RAO) in Socorro, New Mexico, we followed the 
procedures described in IJorstad et al.l (120051 ) to create and analyze the resultant images. We 
refer the polarization position angle measurements to a stable feature in the VLBA polar- 
ized intensity image of the quasar CTA102. We find good agreement between the adjusted 
angles of several objects with archival data and with contemporaneous Very Large Array 
measurements (available on website http://www.vla.nrao.edu/astro/calib/polar/). In our 



analysis, we add 31° to the polarization angle in the pseudocore of 0420—014 to compen sate 
for Faraday rotation, based on the estimate of rotation measure by IJorstad et al.l (120071 ). 



We also observed 0420—014 at two epochs (October 24 and 27) with the Heterodyne 
Spectrometer at the James Clerk Maxwell Telescope (JCMT) at a frequency of 230 GHz. 
The sensitivity of the system was sufficient only to determine the flux density of 0420—014 
relative to the planets Uranus and Mars. 



3. Analysis 



On milliarcsecond (mas) scales, the 43 GHz structure of 0420—014 consists of a compact 
pseudocore and a jet with k nots that propagate away from the pseudocore on paths that 
eventually curve southward (IJorstad et al.ll2005l ). The VLBA images from our campaign 
(Fig. [Q feature a broad jet extending to 0.7 mas south of the pseudocore. The pseudocore, 
with a flux density of ~ 2 Jy, has a low degree of polarization, P < 1%, but a highly variable 
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EVPA (see Fig. [2]). The first distinct jet component, a 0.45 Jy knot lying 0.4 mas south 
of the pseudocore, is more highly polarized (P ~ 7%). Its degree and angle of polarization 
remained steady throughout the campaign to within the uncertainties (cf. Fig. [2]), providing 
a check on the calibration. The very low polarization of the pseudocore at the third epoch 
of VLBA data (Julian date 2453675.9) is an artifact of the long observation time (~ 10 
hours) required to obtain a high-quality image. Upon splitting the observation in two, we 
find that the EVPA varied by ~ 90° from the first half to the second half of the observation, 
thereby cancelling the polarization in the resultant image shown in Figure [TJ The degree of 
polarization of the pseudocore displayed in Figure [2] for this epoch is the average between 
the two halves. For the remaining VLBA images, the EVPA did not change by more than 
2er between the first and second half of the observation. 

The optical polarization underwent pronounced variations in both P and EVPA, as 
shown in Figure [2j Of particular note is that the 43 GHz pseudocore and optical EVPAs 
varied in a similar manner, following each other remarkably well over a range of ~ 80°. 
By eye, we see that a straight line through the three radio points (the polarization of the 
pseudocore was too weak to measure on November 1) provides a reasonable representation 
of the overall optical trend. Upon performing a least-squares linear fit to the EVPA vs. 
time plots that takes into account the uncertainties in the data, we find that the pseudocore 
polarization vector rotates by — 10.5°±0.8° day -1 , while the optical vector rotates by — 11.1°± 
0.3° day -1 . There is no significant time lag between the 43 GHz and optical polarization 
variations. Despite the similarity in EVPA at the two wavebands, the degree of polarization 
P is considerably higher ((P) = 3.0±0.2%) at the optical band than in the 7 mm pseudocore 
((P) = 0.47 ± 0.03%). This implies that ~ 80% of the 43 GHz emission from the pseudocore 
arises from a separate, essentially unpolarized region that is blended with the pseudocore at 
the angular resolution of the VLBA. We suggest that this region is optically thick with a 
tangled magnetic field, in which case the polarization would be very close to zero. A higher 
resolution 86 GHz polarized intensity image should reveal this component. 



4. Discussion 



The striking correspondence between variations in polarization at the two wavebands 
locates the primary site of op tical nonthermal emission in this quasar within the 7 mm 
pseudocore. As emphasized by iMarscherl (J2006bJ), this generally lies well downstream of the 
central engine and "true core," since the self-absorption turnover frequency of the overall 
synch rotron spectrum of 0420—014 occurs at a shorter wavelength, ~ 3 mm (IBloom et al. 
19941 ) . Our flux density measurements of 2.5 Jy at 7 mm and 1.1 Jy at 1.35 mm dur- 



- 5 - 



ing the campaign indeed indicate that the turnover was at a wavelength longer than 1.35 
mm. The true core mig ht be the end of the bulk acceleration and collimation zone (see 
Vlahakis fc Koniglll2004r) beyond which the jet is co nical with essentially constant mean flow 
Lorentz factor (jMarscherlll980l ; iJorstad et al.ll2007l ). 



Magnetic turbulence is the most straightforward way to explain the low fraction al po- 
larization generally observed in the compact jets of quasars at radio wavelengths (e.g., [Jones 
19881 ; iHughesI 120051 ) . We approximate the turbulence as N cells of equal volume, each with 
uniform but randomly oriented magnetic field. The randomness dilutes most of the polar- 
ization of the source such that the mean degree of polarization (P) = P max A^ -1 / 2 , where 
-Pmax = (a + l)/(a + 5/3) is the maximum linear polarization of incoherent synchr otron 
radia t ion and a i s the spectral index of the total flux density spectrum, F u oc u~ a (IBurn 
19661 ; iJoned Il988l ; iHughes fc Miller! 1 19911 ) . The standard deviation about the mean value 
<t(P) = 2~ 1 / 2 (P). The mean fractional optical polarization of ~ 3% then implies that ~ 600 
cells participate in the synchrotron radiation at this waveband. 

The polarization of the pseudocore could vary if there were a changeover in turbulent 
cells involved in the emission during the course of our campaign. This can be accomplished 
i f the radiating electrons are energized by a shock wave that p asses through the turbulence 
( Jones et al.lll985t |Joneall988l ; iMarscher. Gear. &; Travislll992f ) or vice versa. However, only 



in the unlikely case of a face-on propagating shock wave would the partial ordering of the 
magnetic field parallel to the shock front not be apparent from the observer's perspective; 
other wise, the compression would affect the polarization by establishing a preferred EVPA 



(e.g. jHughes fc Miller! 1 1 9 9 lh . Furthermore, electrons accelerated at a shock front lose energy 
as they advect away from the front, with those that attain higher energies — and therefore 
emit at higher frequencies — having shorter radiative lifetimes. This should cause the opti- 
cal emission to occupy a much smaller volume (oc A 1 / 2 ) than does the mm-wave emission 
( IMarscher fc Gearlll985l ). so that the latter should involve ~ 60, 000 cells, or ~ 100 times as 
many as at optical wavelengths. Although this is roughly compatible with the lower mean 
polarization at 7 mm, it cannot explain the similarity in EVPA in the two wavebands, since 
the additional factor of ~ 100 cells should produce a mean polarization angle at 7 mm that 
is essentially independent of the optical EVPA. 

In order for a shock plus turbulence model to be valid, the thickness of the section 
of s hocked plasma imposed by gas dynamics (t he shocked region ends at a rarefaction; see, 
e.g.. lDary fc Marscherlll988l ; lG6mez et al.lll995l ) must be comparable to the radiative lifetime 
of the optically emitting electrons times the speed at which the compressed plasma drifts 
from the acceleration zone at the shock front. The timescale of variability of the optical 
polarization t var « 4 days (see Fig. [2]). If we adopt the apparent speed and angle to 
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the lin e of sight of the jet in 0420—014 derived from VLBA monitoring by iJorstad et al. 



( 120051 ). /9 app ~ 11 and 9 « 3° (Doppler factor of 16 and flow Lorentz factor V « 11), a 
section of jet plasma of thickness x ~ 0.8 pc passes by a stationary shock during a time 
interval of 4 days in the observer's reference frame. The synchrotron radiative lifetime 
of the electrons at optical wavelengths t s ~ 0.18(13/0.1 G) -3 / 2 yr, where we estimate the 
magnetic field B ~ 0.1 G by calculating the value corresp onding to equiparti tion between 
the magnetic and relativistic electron energy densities (see iHoman et al.ll2006l . for evidence 
that the brightness temperatures of pseudocores of blazars are consistent with equipartition). 
The energized electrons travel a distance Tct s ~ 0.6 pc before they can no longer radiate 
in the optical band (the T factor accounts for Lorentz contraction of the jet length in the 
plasma frame). This is essentially the same as the value of x estimated above, so that nearly 
all of the turbulent cells in the emission region are replaced every 4 days, as required by 
the variability. The length of the section compressed by the shock can be less than the 
distance traversed by the electrons radiating at optical wavelengths if the shock results from 
the body pinch-mode instability owing to a n oscillating pressure mi smatch at the boundary 
between the jet and the ex ternal medium (jDalv fc Marscherl Il988l ). As seen in numerical 
gas dynamical simulations (IGomez et al.lll995l ). such shocks have a conical geometry if the 
jet is sufficiently close to being circularly symmetric. The partial ordering of the magnetic 
field owing to compression by a conical sho ck contributes negligibly to the polarization for 
very small viewing a ngles (ICawthornd 120061 ) such as 3°, the value derived for 0420—014 by 
Jorstad et all tood ). 



The ap p arent rotation of the polarization vector can be purely stochastic, as found by 
Jones et al.l (119851 ). who analyzed a model similar to that described above. We have per- 
formed a number of simulations similar to those of Jones et al., incorporating 600 cells, each 
with the same density and strength of a randomly oriented magnetic field. The simulations 
produce artificial plots of percent polarization and EVPA vs. time. Rotations by 100° or 
more occur ~ 10% of the time. Since 0420—014 was only one of 13 objects observed, the 
probability of catching one source at the start of a long rotation episode is reasonably high 
if nearly all blazars behave in a similar manner. 

Our observations could also be explained if a feature, e.g., a shock, with a partially 
ordered magnetic field were to translate down a jet that tw ists by ~ 140° over the lengt h 
of ~ 9 pc that the feature would move in 11 days (see, e.g., iGopal-Krishna fc Wiit all 19921 ). 



This is not implausible, since curved je ts are common in 43 GHz images of blazars (e.g., 



Marscher et al. 



20021 ; I Jorstad et al.ll2005l ). However, such a model does not naturally explain 



the variations in degree of polarization and the fluctuations in EVPA about the rotation trend 
(see Fig. |2|). A similar criticism argues against apparent rotation from changing aberration 
owing to acceleration of a polarized knot as it propagates down the jet (IBlandford &: Konigl 
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19791 ). Furthermore, monthly VLBA monitoring throughout 2006 failed to detect any bright 



knot propagating down the jet that would correspond to a shock wave moving through the 
pseudocore during our campaign. We can also eliminate a simple two-component model in 
which two nearly orthogonally polarized regions vary in relative flux density such that the 
EVPA rotates from that of one feature to that of the other. The maximum rotation possible 
in such a model is less than 90°, while the optical EVPA changed by considerably more than 
100° during the 11 nights of our campaign (see Fig. [2]). 



5. Summary and Conclusions 



Our observation of concurrent polarization position angle at optical wavelengths and 7 
mm as it swings by ~ 80° over an interval of nine days demonstrates that the main region 
producing polarized optical emission is the same as that giving rise to the pseudocore at 43 
GHz. During this rotation, the degree of polarization varied and the EVPA fluctuated about 
the mean trend in a manner consistent with a chaotic magnetic field. We propose a model 
that explains the variations via changeovers in the ~ 600 turbulent cells involved in the 
emission over a period of a few days. This is accomplished as the plasma advects through a 
standing, conical shock where electrons are accelerated at the front. The polarized emission 
at 43 GHz must originate in the same volume as does the optical radiation, a situation that 
requires the region compressed by the shock to be limited by gas dynamics rather than by 
radiative energy losses. This agrees with gas dynamical simulations th at show alternatin g 
regions of conical shocks and rarefactions in jets with circular symmetry (IGomez et al.lll995l ). 
We therefore conclude that the 43 GHz pseudocore in 0420—014 corresponds to a standing 
shock system rather than the site where the optical depth is roughly unity. This is likely to 
be the case in m any blazars, since st ationary features near the pseudocore in parsec-scale 
jets are common (IJorstad et al.ll200ll ). 



We are in the process of analyzing data from other objects observed during the campaign 
as well as objects from a similar campaign in 2006 March-April. These data will determine 
how frequently rapid, correlated variations in optical and high-frequency radio EVPAs occur 
in blazars. If this behavior is common, then we can search for similarity in polarization in a 
number of blazars to identify the locations of optical synchrotron emission on VLBI images. 
This would allow us to use correlations between optical or mm-wave and X-ray as well as 
7-ray variations in total flux density to determine the location ( s) of nontherma l emission 
across a factor of ~ 10 12 in frequency (IMarscher fc Jorstadll2005l ; iMarscherl l2006al ) . 



This material is based on work at Boston University supported by the National Science 



- 8- 



Foundation under grant AST-0406865, and at St. Petersburg State University supported 
by grant no. 05-02-17562 from the Russian Fund for Basic Research. P.S.S. acknowledges 
support from NASA contract 1256424. The VLBA is an instrument of the National Radio 
Astronomy Observatory, a facility of the National Science Foundation operated under co- 
operative agreement by Associated Universities, Inc. The JCMT is operated by the Joint 
Astronomy Centre on behalf of the Particle Physics and Astronomy Research Council of the 
United Kingdom, the Netherlands Organization for Scientific Research, and the National 
Research Council of Canada. 

REFERENCES 

Blandford, R. D., & Konigl, A. 1979, ApJ, 232, 34 
Bloom, S. D., et al. 1994, AJ, 108, 398 
Burn, B. J. 1966, MNRAS, 133, 67 
Cawthorne, T. V. 2006, MNRAS, 367, 851 
Cawthorne, T. V., & Cobb, W. K. 1990, ApJ, 350, 536 
Daly, R. A., & Marscher, A. P. 1988, ApJ, 334, 539 

Gabuzda, D. C, Rastorgueva, E. A., Smith, P. S., & O'Sullivan, S. P. 2006, MNRAS, 369, 
1596 

Gabuzda, D. C, & Sitko, M. L. 1994, AJ, 107, 884 

Gomez, J. L., Marti, J. M., Marscher, A. P., Ibahez, J. M., & Marcaide, J. M. 1995, ApJ, 
449, L19 

Gopal-Krishna, & Wiita, P. J. 1992, A&A, 259, 109 
Hewitt, A., & Burbidge, G. 1993, ApJS, 87, 451 
Homan, D. C, et al. 2006, ApJ, 642, LI 15 

Hughes, P. A., & Miller, L. 1991, in Beams and Jets in Astrophysics, ed. P. A. Hughes, 
Cambridge Astrophysics Series, (Cambridge Univ. Press), 1 

Hughes, P. A. 2005, ApJ, 621, 635 

Jones, T. W. 1988, ApJ, 332, 678 



- 9- 



Jones, T. W., et al. 1985, ApJ, 290, 627 
Jorstad, S. G., et al. 2001, ApJS, 134, 181 
Jorstad, S. G., et al. 2005, AJ, 130, 1418 
Jorstad, S. G., et al. 2007, A J, submitted 
Kdnigl, A. 1981, ApJ, 243, 700 
Lister, M. L., & Smith, P. S. 2000, ApJ, 541, 66 
Marscher, A. P. 1980, ApJ, 235, 386 

Marscher, A.P, Jorstad, S.G., Mattox, J.R., & Wehrle, A.E. 2002, ApJ, 577, 85 

Marscher, A. P. 2006a, in Blazar Variability Workshop II: Entering the GLAST Era, ed. H. 
R. Miller et al., ASP Conf. Ser., 350, 155 

Marscher, A. P. 2006b, in Relativistic Jets: The Common Physics of AGN, Microquasars 
and Gamma-Ray Bursts, ed. P.A. Hughes & J.N. Bregman, AIP Conf. Ser., 856, 1 

Marscher, A. P., & Gear, W. K. 1985, ApJ, 298, 114 

Marscher, A. P., Gear, W. K., & Travis, J. P. 1992, in Variability of Blazars, ed. E. Valtaoja 
& M. Valtonen, Cambridge U. Press, 85 

Marscher, A. P., & Jorstad, S. G. 2005, in Astronomical Polarimetry: Current Status and 
Future Directions, ed. A. Adamson et al., ASP Conf. Ser., 343, 469 

Schmidt, G. D., Stockman, H. S., & Smith, P. S. 1992, ApJ, 398, L57 

Sitko, M. L., Schmidt, G. D., & Stein, W. A. 1985, ApJS, 59, 323 

Smith, P. S., & Balonek, T. J. 1998, PASP, 110, 1164 

Smith, P. S., Schmidt, G. D., Hines, D. C, & Foltz, C. B. 2003, ApJ, 593, 676 

Vlahakis, N., & Kdnigl, A. 2004, ApJ, 605, 656 

Wardle, J. F. C, & Kronberg, P. P. 1974, ApJ, 194, 249 

Wardle, J. F. C, Moore, R. L., & Angel, J. R. P. 1984, ApJ, 279, 93 



This preprint was prepared with the AAS IATgX macros v5.2. 



-10- 



24 Oct. 05 28 Oct. 05 1 Nov. 05 2 Nov. 05 




1 2 3 4 5 



Fig. 1. — 43 GHz (7 mm) images of 0420-014 at the four epochs of VLBA observations. 
We refer to the region surrounding the peak in total intensity as the "pseudocore." The 
sacle of the horizontal, right ascension axis, are the same as for the vertical, declination 
axis. Contours represent 1, 2, 4, 8, 16, 32, 64, and 90% of the peak total intensity of 2.0 
Jy/beam, while the sticks indicate polarized intensity (length of 0.1 mas corresponds to 6.2 
mJy/beam) and direction of electric vector. Correction for statistical bias and for Faraday 
rotation toward the pseudocore of —31° at 43 GHz has not been applied to the displayed 
polarization vectors. As noted in the text, the very low polarization of the pseudocore in 
the 2005 November 1 image is an artifact of a pronounced change in polarization position 
angle during the ten-hour observation. The rms polarization noise level is 1.8 mJy/beam. 
Dimensions of the restoring beam (shown as a shaded ellipse in the lower right corner) are 
0.40x0.17 mas along position angle —11°. 
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Fig. 2. — Top: Variation of percent polarization with time of 0420—014 in the optical band 
(450-750 nm; filled circles and solid line), in the 7 mm pseudocore (open diamonds and 
dotted line), and in the first jet component in the 7 mm image (x and dashed line). The 
degrees of polarization are corrected for statistical bias (IWardle fc Kronberglll974l ). Bottom: 
Variation of EVPA with time of the same three components of PKS 0420—014. The EVPA 
is not plotted when P/a(P) < 1.5 or for the pseudocore at JD 2453675.9 when it varied 
greatly during the VLBA observation. A correction of 31° has been added to the EVPAs 
of the pseudocore to compensate for Faraday rotation. For reference, JD 2453670 is 2005 
October 26. 



